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Abstract 

A cosmic acceleration mechanism is introduced which is based on the wakefields 
excited by the Alfven shocks in a relativistically flowing plasma, where the energy 
gain per distance of a test particle is Lorentz invariant. We show that there exists 
a threshold condition for transparency below which the accelerating particle is 
collision-free and suffers little energy loss in the plasma medium. The stochastic 
encounters of the random accelerating-decelerating phases results in a power- 
law energy spectrum: /(e) oc 1/e 2 . The environment suitable for such plasma 
wakefield acceleration can be cosmically abundant. As an example, we discuss 
the possible production of super-GZK ultra high energy cosmic rays (UHECR) 
through this mechanism in the atmosphere of gamma ray bursts. We show that 
the acceleration gradient can be as high as G ~ 10 16 eV/cm. The estimated event 
rate in our model agrees with that from UHECR observations. 
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Ultra high energy cosmic ray (UHECR) events exceeding the Greisen-Zatsepin-Kuzmin 
(GZK) cutoff[l] (5 x 10 19 eV for protons originated from a distance larger than ~ 50 Mps) 
have been found in recent years[2, 3, 4, 5]. Observations also indicate a change of the power- 
law index in the UHECR spectrum (events/area/time oc e~ a ) from a ~ 3 to a smaller value, 
at energy around 10 18 — 10 19 eV (the so-called "ankle"). These present an acute theoretical 
challenge regarding their composition as well as their origin [6]. 

So far the theories that attempt to explain the UHECR can be largely categorized into 
the "top-down" and the "bottom-up" scenarios. In addition to relying on exotic particle 
physics beyond the standard model, the main challenges of top-down scenarios are their 
difficulty in compliance with the observed event rates and the energy spectrum[6], and the 
fine-tuning of particle lifetimes. The main challenges of the bottom-up scenarios, on the other 
hand, are the GZK cutoff, as well as the lack of an efficient acceleration mechanism [6] . To 
circumvent the GZK limit, several authors propose[7] that it is neutrinos, instead of protons, 
that are the actual messenger. With the much weaker (electroweak) interaction, neutrinos 
can propagate across a much larger cosmic distance than the GZK radius. These ultra-high 
energy neutrinos would eventually be annihilated by cosmic background neutrinos and turn 
into Z-bosons. If this happens in our surrounding cluster, the proton that is produced from 
the Z-decay can reach the Earth. For such a cascade scenario to work, it requires that the 
original particle, say protons, be several orders of magnitude more energetic than the one 
eventually reaches the Earth. 

Even if the GZK-limit can be circumvented through the Z-burst, the challenge for a vi- 
able acceleration mechanism remains, or becomes even more acute. This is mainly because 
the existing paradigm for cosmic acceleration, namely the Fermi mechanism [8] , as well as its 
variants, such as the diffusive shock acceleration [9], are not effective in reaching ultra high 
energies[10]. These acceleration mechanisms rely on the random collisions of the high energy 
particle against magnetic field domains or the shock media, which necessarily induce increas- 
ingly more severe energy losses at higher particle energies. There is therefore a "diminishing 
return" where the energy gain becomes evermore difficult, if not impossible. 

From the experience of terrestrial particle accelerators, we learn that it takes several 
qualifications for an accelerator to operate effectively. First, the particle should gain energy 
through the interaction with the longitudinal electric field of a subluminous (v < c) electro- 
magnetic (EM) wave. In such a setting the accelerated particle can gain energy from the 
field over a macroscopic distance, much like how a surfer gains momentum from an ocean 
wave. It is important to note that such a longitudinal field is Lorentz invariant, meaning 
that the acceleration gradient is independent of the instantaneous energy of the accelerating 
particle. Second, such a particle-field interaction should be a non-collisional process. This 
would help to avoid severe energy loss through inelastic scatterings. Third, to avoid ex- 
cessive synchrotron radiation loss, which scales as particle energy squared, the accelerating 
particle should avoid any drastic bending beyond certain energy regime. We believe that 
these qualifications for terrestrial accelerators are also applicable to celestial ones. 

Although they are still in the experimental stage, the "plasma wakefield accelerator" 
concepts[ll, 12], promise to provide all the conditions stated above. Plasmas are capable of 
supporting large amplitude electro-static waves with phase velocities near the speed of light. 
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Such collective waves, or "wakefields" , can be excited by highly concentrated, relativistic EM 
energies such as lasers[ll] and particle beams[12]. A trailing particle can then gain energy 
by riding on this wakefield. Although hard scatterings between the accelerating particle 
and the plasma medium is inevitable, under appropriate conditions, as we will demonstrate 
below, the particle can be collision-free. 

In this Letter we demonstrate that magneto-shocks (Alfven shocks) in a relativistic 
plasma flow can also excite large amplitude plasma wakefields, which in turn can be highly 
efficient in accelerating ultra high energy particles. We note that the physical conditions re- 
quired for such acceleration are cosmically abundant and therefore generic. When applying 
this mechanism to the atmosphere of gamma ray bursts (GRBs), we show that protons can 
be accelerated to energies much beyond ZeV (10 21 eV) with a power-law spectrum. 

It is well-known that an ordinary Alfven wave propagating in a stationary magnetized 
plasma has a velocity 

which is typically much less than the speed of light. Here Bo is the longitudinal magnetic 
field and n p is the density of the magnetized plasma. Such a slow wave is ordinarily not 
suitable for accelerating relativistic particles. The situation changes when the plasma as a 
whole moves with a relativistic bulk velocity V p < c. The standard method of obtaining the 
linear dispersion relation of waves in a magnetized plasma leads to 

T p (u - V p k ± u Bi /T p ) (u - V p k T u;Be/T p ) ' 1 ' 
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where k and uo are the wave number and the frequency of the EM wave, respectively, io pi)Pe = 
(A-Ke 2 n p /mi :e ) 1 ^ 2 are the plasma frequencies for ions and electrons, and ujBi,Be = (e-Bo/^i.e) 1 ^ 2 
are the ion and electron cyclotron frequencies. Here T p is the Lorentz factor of the bulk 
plasma flow. Figure 1 shows the dispersion relations of various transverse EM waves that 
propagate along the direction of Bo with and without the plasma bulk flow V p . In Fig. 
1(a) we see that outside the lightcone (superluminous, or v p h > c) lie the regular EM 
waves, whose asymptotic dispersion is u = kc. Within the lightcone (subluminous), there 
are two additional branches, the whistler wave (an electron branch mode) and the Alfven 
wave whose frequency remains quite low and its electric field is much smaller than the 
magnetic one, i.e., E/B = v A /c <C 1 in the absence of flow. In the case where the bulk 
flow of the plasma approaches the speed of light, however, both the whistler and Alfven 
waves now acquire a phase velocity close to c from below (see Fig. 1(b)) and enhances 
the ratio of E/B to ~ V p /c < 1. It is important to emphasize that as the Alfven wave 
dispersion approaches the lightcone, it becomes indistinguishable from a bona fide EM wave, 
and is therefore able to excite relativistic plasma wakefields, much like lasers and particle 
beams. Preliminary results from simulations indicate that such relativistic Alfven waves can 
indeed excite plasma wakefields [13]. It is also important to realize that in this relativistic 
flow the excited wakefields are all in one direction, which contributes to the unidirectional 
acceleration. With our applications to astrophysical problems in mind, the Alfven-wave- 
plasma interaction relevant to us is in the nonlinear regime. 
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The plasma wakefield in the nonlinear regime has been well-studied[14]. The nonlinearity 
is determined by the driving EM wave's ponderomotive potential, which is governed by its 
normalized vector potential ao = eE/mcu. When this parameter exceeds unity, nonlinearity 
is strong [11] so that additional important physics incurs. For a stationary plasma, the 
maximum field amplitude that the plasma wakefield can support is 

m e cujp 

E max ~ -C/wbOo = o-o , (3) 

e 

which is enhanced by a factor ao from the cold wavebreaking limit (the naively assumed 
maximum field), _E w b = m e cujp/e, of the linear regime. In a relativistic plasma flow with 
a Lorentz factor T p , the cold wavebreaking field is reduced by a factor T^ 2 due to Lorentz 
contraction. The maximum "acceleration gradient" G experienced by a singly-charge particle 
riding on this plasma wakefield is then 



G = eE' max ^a m e c^^^. (4) 

The plasma wavelength, in the mean time, is stretched also by a factor ao from that in the 
linear regime. So in a plasma flow the wavelength is 




ApAr — —OoK, ~ a W , (5) 



where r e = e 2 /m e c 2 = 2.8 x 10~ 13 cm is the classical electron radius. 

To determine the maximum possible energy gain, we need to know how far can a test 
particle be accelerated. At ultra high energies once the test particle encounters a hard scat- 
tering or bending, the hard-earned kinetic energy would most likely be lost. The scattering of 
an ultra high energy proton with the background plasma is dominated by the proton-proton 
collision. Existing laboratory measurements of the total pp cross section scales roughly as 
c PP = o~o • {1 + 6.30 x 10~ 3 [log(s)] 2 - 1 }, where a 32mb and the center-of-mass energy- 
squared, s, is given in (GeV) 2 . In our system, even though the UHE protons are in the ZeV 
regime, the center-of-mass energy of such a proton colliding with a comoving background 
plasma proton is in the TeV range, so it is safe to ignore the logarithmic dependence and 
assume a constant total cross section, a pp ~ er ~ 30 mb in the ZeV energy regime. Since in 
astrophysical settings an out-bursting relativistic plasma dilutes as it expands radially, its 
density scales as n p (r) = n p0 (Ro/r) 2 , where n p0 is the plasma density at a reference radius 
Rq . The proton mean-free-path can be determined by integrating the collision probability 
up to unity, 

rRo+L mfp o-ppWpCr) _ rRo+L mip a pp np R 2 

' Rq J Rq .i J ; 

We find 



= Z-^+^p <Tppnp( r ) dr = f Ro+L ^ ytgo R 2 dr ^ 

jRn r„ jRn r„ r 2 



CTppUpoRi 



Rn 



Tp L Ro + L m f p 



(7) 
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Since L mfp is positive definite, < [1 — R /(R + L m{p )] < 1. Therefore the solution to 
L m f p does not exist unless the coefficient, a pp n P QR /T p > 1. That is there exists a threshold 
condition below which the system is collision-free: 

(T pp n p oRo = j ^ 

When a system is below this threshold, a test particle can in principle be accelerated un- 
bound. In practice, of course, other secondary physical effects would eventually intervene. 

In a terrestrial accelerator, the wakefields are coherently excited by the driving beam, 
and the accelerating particle would ride on the same wave crest over a macroscopic distance. 
There the aim is to produce near-monoenergetic final energies (and tight phase-space) for 
high energy physics and other applications. In astrophysical settings, however, the drivers, 
such as the Alfven shocks, will not be so organized. A test particle would then face random 
encounters of accelerating and decelerating phases of the plasma wakefields excited by Alfven 
shocks. 

The stochastic process of the random acceleration-deceleration can be described by the 
distribution function f(e,t) governed by the Chapman-Kolmogorov equation[15, 16] 

f) r+oo r+oo 

-f= / d(Ae)W(e-Ae,Ae)f(e-Ae,t)- / d(Ae)W(e, Ae)/(e, t) - u(e)f(e, t) . (9) 

OT J—oo J — oo 

The first term governs the probability per unit time of a particle "sinking" into energy e 
from an initial energy e — Ae while the second term that "leaking" out from e. The last 
term governs the dissipation due to collision or radiation, or both. As we will demonstrate 
later, the astrophysical environment that we invoke for the production of UHECR is below 
the collision threshold condition, and so accelerating particles are essentially collision-free. 

The radiation loss in our system is also negligible. As discussed earlier, in a relativis- 
ts flow the transverse E and B fields associated with the Alfven shock are near equal in 
magnitude. Analogous to that in an ordinary EM wave, an ultra relativistic particle (with 
a Lorentz factor 7) co-moving with such a wave will experience a much suppressed bending 
field, by a factor I/7 2 . Furthermore, the plasma wakefield acceleration takes place in the 
region that trails behind the shock (and not in the bulk of the shock) where the acceler- 
ating particle in effect sees only the longitudinal electrostatic field colinear to the particle 
motion[14]. We are therefore safe to ignore the radiation loss entirely as well. We can thus 
ignore the dissipation term in the Chapman-Komogorov equation and focus only on the 
purely random plasma wakefield acceleration-deceleration. 

Assuming that the energy gain per phase encounter is much less than the final energy, 
i.e., Ae <C e, we Taylor-expand W(e — Ae,Ae)f(e — Ae) around W(e, Ae)/(e) in the sink 
term and reduce Eq.(9) to the Fokker-Planck equation 



d d r +OD d 2 r+°° Ae 2 

dt f = d~6 J 00 ^ Ae ^ (C ' ^ + J 00 d ^) — W ^ Ae )^ 6 ' *) • ( 10 ) 



We now assume the following properties of the transition rate W(e, Ae) for a purely 
stochastic process: 
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a) W is an even function; 

b) W is independent of e; 

c) W is independent of Ae. 

Property a) follows from the fact that in a plasma wave there is an equal probability of 
gaining and losing energy. In addition, since the wakefield amplitude is Lorentz invariant, 
the chance of gaining a given amount of energy, Ae, is independent of the particle energy e. 
Finally, under a purely stochastic white noise, the chance of gaining or losing any amount 
of energy is the same. Based on these arguments we deduce that 

where r is the typical time of interaction between the test particle and the random waves 
and G is the maximum acceleration gradient (cf. Eq.(4)). We note that there is a stark 
departure of the functional dependence of W in our theory from that in Fermi's mechanism, 
in which the energy gain Ae per encounter scales linearly and quadratically in e for the 
first-order and second-order Fermi mechanism, respectively. 

To look for a stationary distribution, we put df/dt = 0. Since W is an even function, the 
first term on the RHS in Eq.(10) vanishes. To ensure the positivity of particle energies before 
and after each encounter, the integration limits are reduced from (—00, +00) to [— e, +e], and 
we have 

d 2 f+£ Ae 2 

^ J e d(Ae)—W(e,Ae)f(e)=Q. (12) 

Inserting W from Eq.(ll), we arrive at the energy distribution function that follows power- 
law scaling, 

/(e) = , (13) 

where the normalization factor eo is taken to be the mean energy of the background plasma 
proton, eo ~ Y p m p c 2 . The actually observed UHECR spectrum is expected to be degraded 
somewhat from the above idealized, theoretical power-law index, a — 2, not only due to pos- 
sible departure of the reality from the idealized model, but also due to additional intermediate 
cascade processes that transcend the original UHE protons to the observed UHECRs. 

We note that a power-law energy spectrum is generic to all purely stochastic, collisionless 
acceleration processes. This is why both the first and the second order Fermi mechanisms 
also predict power-law spectrum, if the energy losses, e.g., through inelastic scattering and 
radiation (which are severe at ultra high energies), are ignored. The difference is that 
in the Fermi mechanism the stochasticity is due to random collisions of the test particle 
against magnetic walls or the shock medium, which necessarily induce reorientation of the 
momentum vector of the test particle after every diffusive encounter, and therefore should 
trigger inevitable radiation loss at high energies. The stochasticity in our mechanism is due 
instead to the random encounters of the test particle with different accelerating-decelerating 
phases. As we mentioned earlier, the phase vector of the wakefields created by the Alfven 
shocks in the relativistic flow is nearly unidirectional. The particle's momentum vector, 
therefore, never changes its direction but only magnitude, and is therefore radiation free in 
the energy regime that we consider for proton acceleration. 
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We now apply our acceleration mechanism to the problem of UHECR. GRBs are by far 
the most violent release of energy in the universe, second only to the big bang itself. Within 
seconds (for short bursts) about cgrb ~ 10 52 erg of energy is released through gamma rays 
with a spectrum that peaks around several hundred keV. Existing models for GRB, such 
as the relativistic fireball model [17], typically assume neutron-star-neutron-star (NS-NS) 
coalescence as the progenitor. Neutron stars are known to be compact (i?NS ~ O(lO)km) 
and carrying intense surface magnetic fields (B NS ~ 10 12 G). Several generic properties are 
assumed when such compact objects collide. First, the collision creates sequence of strong 
magneto-shocks (Alfven shocks). Second, the tremendous release of energy creates a highly 
relativistic out-bursting fireball, most likely in the form of a plasma. 

The fact that the GRB prompt (photon) signals arrive within a brief time-window implies 
that there must exists a threshold condition in the GRB atmosphere where the plasma 
becomes optically transparent beyond some radius i?o from the NS-NS epicenter. Applying 
Eq. (8) to the case of out-bursting GRB photons, this condition means 



where a c = (7rr 2 )(m e /o;GRB)[log(2a;GRB/''^e) + l/2] ~ 2x 10~ 25 cm 2 is the Compton scattering 
cross section. Since o pp < a c , the UHECRs are also collision-free in the same environment. 
There is clearly a large parameter space where this condition is satisfied. To narrow down 
our further discussion, it is not unreasonable to assume that Ro ~ O(10 4 )km. A set of self- 
consistent parameters can then be chosen: n p0 ~ 10 20 cm~ 3 , T p ~ 10 4 , and e ~ 10 13 eV = ei 3 . 

To estimate the plasma wakefield acceleration gradient, we first derive the value for the 
cto parameter. We believe that the megneto-shocks constitute a substantial fraction, say 
i] a ~ 10~ 2 , of the total energy released from the GRB progenitor. The energy Alfven shocks 
carry is therefore 6a ~ 10 50 erg. Due to the pressure gradient along the radial direction, 
the magnetic fields in Alfven shocks that propagate outward from the epicenter will develop 
sharp discontinuities and be compactified[18]. The estimated shock thickness is ~ 0(l)m at 
i?o ~ O(10 4 )km. From this and one can deduce the magnetic field strength in the Alfven 
shocks at Rq, which gives Ba ~ 10 10 G. This leads to ao = eEU/mca^ ~ 10 9 . Under these 
assumptions, the acceleration gradient G (cf. Eq.(4)) is as large as 



Although the UHE protons can in principle be accelerated unbound in our system, the 
ultimate maximum reachable energy is determined by the conservation of energy and our 
assumption on the population of UHE protons. Since it is known that the coupling between 
the ponderomotive potential of the EM wave and the plasma wakefield is efficient, we assume 
that the Alfven shock energy is entirely loaded to the plasma wakefields after propagating 
through the plasma. Furthermore, we assume that the energy in the plasma wakefield is 
entirely reloaded to the UHE protons through the stochastic process. Thus the highest 
possible UHE proton energy can be determined by energy conservation 




(14) 




(15) 




(16) 
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which gives 

e m = e 13 exp(r] a e GRB /N VHE e 13 ) . (17) 

This provides a relationship between the maximum possible energy, e m , and the UHE proton 
population, A^he- We assume that r] b ~ 1CT 2 of the GRB energy is consumed to create the 
bulk plasma flow, i.e., t^cgrb ~ N p Y p m p c 2 ~ N p e 13 , where N p is the total number of plasma 
protons. We further assume that r] c ~ 10~ 2 of the plasma protons are trapped and accelerated 
to UHE, i.e., iVuHE ~ VcN p . Then we find e m ~ ei 3 exp^/r^T^). We note that this estimate 
of e m is exponentially sensitive to the ratio of several efficiencies, and therefore should be 
handled with caution. If the values are indeed as we have assumed, r] a /r] b r] c ~ O(10 2 ), then 
e m is effectively unbound until additional limiting physics enters. Whereas if the ratio is 
~ O(10) instead, the UHE cannot even reach the ZeV regime. The validity of our assumed 
GRB efficiencies then relies on the consistency check against observations. 

In addition to the energy production issue, equally important to a viable UHECR model 
is the theoretical estimate of the UHECR event rates. The NS-NS coalescence rate is believed 
to be about 10 events per day in the entire Universe[19, 20]. This frequency is consistent 
with the observed GRB events, which is on the order of /grb ~ 10 35 per year. 

In the Z-burst scenario an initial neutrino energy above 10 21 eV[7] or 10 23 eV[21] is required 
(depending on the assumption of the neutrino mass) to reach the Z-boson threshold. For 
the sake of discussion, we shall take the necessary neutrino energy as e v > 10 22 eV. Such 
ultra high energy neutrinos can in principle be produced through the collisions of UHE 
protons with the GRB background protons: pp^ir + X^fi + v + X. All UHE protons 
with energy e>22 > 10 22 eV should be able to produce such neutrinos. The mean energy 
(by integrating over the distribution function /(e)) of these protons is (e >2 2) ~ O(100)e22- 
Therefore the multiplicity of neutrinos per UHE proton is around H{ p ^ u ) ~ O(10) — O(100). 
At the opposite end of the cosmic process, we also expect multiple hadrons produced in 
a Z-burst. The average number of protons that Z-boson produces is ~ 2. 7 [22]. Finally, 
the population of UHE protons above 10 22 eV is related to the total UHE population by 

A^>22 ~ (£13/^22) AuHE ~ 

Putting the above arguments together, we arrive at our theoretical estimate of the ex- 
pected UHECR event rate on earth, 

1 

Auhecr(> 10 20 eV) = fGRBH(p->u)H(z->p)N >2 2-r-^2 — 

47r - K GRB 

~ fGRBV(p^v)V(Z^p)VbVc €GRB -j p2 • (18) 

e 2 2 4:7lK GRB 

The typical observed GRB events is at a redshift z ~ 0(1), or a distance -Rgrb ~ 10 22 km. 
Our estimate of observable UHECR event rate is therefore 

AWm(> 10 20 eV) = O(l)/100km 2 /yr/sr , (19) 

which is consistent with observations, or in turn this observed event rate can serve as a 
constraint on the various assumptions of our specific GRB model. 

We have demonstrated that plasma wakefields excited by Alfven shocks in a relativistic 
plasma flow can be a very efficient mechanism for cosmic acceleration, with a power-law 
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energy spectrum. When invoking GRBs as the sites for UHECR production with a set of 
reasonable assumptions, we show that our estimated UHECR event rate is consistent with 
observations. This cosmic acceleration mechanism is generic, and can in principle be applied 
to other astrophysical phenomena, such as blazars[23\. It is generally believed that the AGN 
jets are relativistic plasmas. The observed "lumps", or density concentrations, in the jet 
may well serve as the driver to excite plasma wakefields. These wakefields can accelerate 
electrons as well as protons to multi-TeV energies. Bent by the confining helical magnetic 
fields in the jet, these high energy electrons can radiate hard photons in the TeV range, 
while the protons can cascade into high energy neutrinos. We will present a more detailed 
discussion on blazars in a separate paper. 

We appreciate helpful discussions with J. Arons, R. Blandford, P. Meszaros. 
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Figure 1: The dispersion relations for stationary and relativistic plasma flows. 
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